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ABSTRACT 



The classic single echelon inventory model is restruc- 
tured as a two echelon problem in which demand and resupply 
are deterministic. Using cost minimization as the objective 
three models are developed which address the problems of 
(1) no stockouts allowed (EOQ) , (2) backorders allowed, and 

(3) finite production with no stockouts allowed. General 
solutions for the optimal policy are obtained in the EOQ and 
finite production models. In the backorder model, the ana- 
lytical argument is limited to the case in which only time 
dependent backorder costs occur. Algorithms are developed 
for solving problems for all three models, and selected para 
meter values are used to test the behavior of the models. 
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I . 



INTRODUCTION 



Traditional inventory analysis addresses the single 
echelon inventory problem. In this problem, customer demands 
are received and issues are made from a single outlet. The 
outlet in turn is replenished from a single source. 

Unfortunately, most real world inventory systems are not 
this simple. A major manufacturer, for example, usually has 
a production and distribution system that includes (1) fac- 
tories, (2) factory warehouses, (3) regional warehouses, 

(4) local warehouses, and (5) retail outlets [Hadley and 
Whitin] . Such a system is an example of a multiechelon inven- 
tory system. Each level of the system is a separate, distinct 
echelon. Figure 1 depicts graphically this organization. The 
United States Navy with its organization of inventory control 
points, supply centers, supply depots, shop stores, tenders and 
underway replenishment ships provides an example of a large 
multiechelon inventory system. 

In defending their single echelon models Hadley and Whitin 
argue that, even though most real world systems are multi- 
echelon, it is often true that the system need not be treated 
as multiechelon. They contend that a different organization 
frequently operates each level (echelon) of the system. In 
the example of Figure 1, a factory and its associated ware- 
houses might be operated by an equipment manufacturer whereas 
the regional warehouses might be separate privately owned 
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Figure 1. Typical Multiechelon Inventory System 



distributorships not under the administrative control of the 
manufacturer. In like manner, the local outlets might be 
individually owned retail stores. Thus we would have a sys- 
tem comprised of several single echelons which are linked by 
a physical dependence but separated by administrative 
controls . 

While it is clear that there are many examples of multi- 
echelon systems that can be treated as a series of independent 
single echelon problems, it should be equally clear that there 
is a large class of systems which must be treated as strongly 
dependent multiechelon. 

It is the objective of this thesis to develop and investi- 
gate the behavior of multiechelon models constructed from the 
classic single echelon inventory models. Specifically, the 
following deterministic single-echelon models of Hadley and 
Whitin will be restructured as two-echelon models: 
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a. No stockouts allowed (Economic Order Quantity) , 

b. Backorders allowed, 

c. Finite production with no stockouts allowed. 

Like their single echelon counterparts, these two-echelon 
models will be single item, single source models. 
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II. THE TWO ECHELON MODEL WITH NO STOCKOUTS 



In any study of inventory theory, one of the first models 
investigated is the deterministic demand model with no stock- 
outs allowed; often referred to as the Economic Order Quantity 
(EOQ) model. While it is true that complete deterministic 
demand is almost never known, it is felt that the mathematics 
of this model provide a good starting point for any inventory 
analysis. Further, it is felt that the deterministic model 
will provide an insight into the operation of a stochastic 
demand model. It is considered appropriate, therefore, to 
begin this multiechelon analysis by considering a two echelon 
extension of the classic EOQ model. 

A. MODEL FORMULATION 

The two echelon model is based on the following 
assumptions : 

1. The upper echelon always replenishes its supply from 
the same outside source. The lower echelon always replenishes 
its supply from the upper echelon. 

2. The upper echelon is always able to meet the demand 
of the lower echelon. 

3. External customer demand always occurs at the lower 
echelon. 

4. The external demand rate is deterministic, continuous, 
and constant with a value of A units per year. 
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5. Neither backorders nor lost sales are allowed. 

6. Procurement lead time (PLT) is negligible. 

7. The order quantity at the upper echelon, Q2, is an 
integer multiple of the lower echelon order quantity, Ql. 

Let n = (^/Qj ; then this assumption requires n > l and 
integer. 

The optimal policy throughout this thesis will be that 
which minimizes average annual variable system cost, subject 
to the constraint (s) of the model. The form of the average 
annual variable system cost will be developed by first deter- 

« 

mining the system cost per cycle, and then dividing this cost 
by the cycle length. The total system variable cost per cycle 
is the sum of the individual echelon variable cycle costs, 
K^(c) and (c) . The cycle length is defined as the time be- 
tween receipt of two successive orders at the upper echelon 
(echelon two) . 

The cost per cycle at echelon two is comprised of an 
ordering cost, an inventory holding cost (IHC) and a purchase 
cost. The ordering cost is assumed to be a constant cost per 
order which includes the administrative costs associated with 
inventory review and order (contract) preparation. This cost 
is independent of the quantity on hand or on order. The IHC 
includes a warehousing cost, an obsolescence cost and a fore- 
gone opportunity cost. It is assumed to be a function of the 
inventory on hand and to be expressable by 

T 

IHC = I C fci 9 (t)dt, 

0 
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where is the average cost per dollar invested in inventory 
per unit time, C is the item unit cost, T is the cycle length, 
and ( t ) is the on hand inventory at echelon two at time t. 

The purchase cost is assumed to be independent of the quantity 
produced, and can be expressed as C^C where Q 2 is the quantity 
procured . 

The variable cost per cycle at echelon two can now be 
expressed as 
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( 2 - 1 ) 



K 2 (c) 



T 

A 2 + CQ 2 + I 2 C 2 ( t) dt 



The height of the on-hand inventory at echelon two is 
Q 2 - Q 1 = [n-1] Q ± 



at the beginning of the cycle since Q 2 is ordered, arrives 
immediately, and is partially used to meet an immediate demand 
of Q 1 units. The next demand at echelon two is for an amount 
which occurs Q^/A time units later. Similarly, another 
demand of occurs after Q^/A additional time units have 
passed. The on-hand level at echelon two continues to de- 
crease in steps of until it reaches zero. It remains at 
zero for the last Q-^/x time units of T. The area under the 
on-hand inventory curve for one cycle of echelon two is 
therefore 





n-1 

l mQ 

m=l 



1 



n [n-1] ^2 
2A y l * 



From (2-1) , the total variable cost per cycle at echelon 
two is 



K 2 (c) - A 2 + CQ 2 + 



I 2 C n [n-1] 



2A 



( 2 - 2 ) 



From the assumption that Q 2 = nQ 1 , where n is integer, 
it follows that the cycle length at echelon two is equal to n 
reorders- at echelon one. Therefore, from Figure 2, 



nQ, 



T = 



(2-3) 



The average annual cost at echelon two is obtained by 
dividing (2-2) by (2-3) which yields 
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K. 



A 9 A I C[n-1]Q, 

^ + CA + • 



2Q, 



or, as a function of Q-^ and n, 



K, 



A„A I CQ, [n-1] 

+ CA + — = 

nQ, 2 



(2-4) 



Development of the cost per cycle at echelon one follows 
an argument that is analogous to the echelon two development, 
except that there is no purchase cost incurred at echelon one. 
Also, from the assumption that Q 2 is an integer multiple of 
, it follows that in each cycle there are n reorders at 
echelon one, i.e., K^(c) is linear in n. Thus the cost per 
cycle can be written as 

nl 1 CQ? 

K^c) = nh 1 + — ~ — — . (2-5) 



When (2-5) is divided by (2-3) the average annual variable 
cost is 

A.A I,CQ, 

K i = 07 + Hr- 1 • (2 " 6) 

The total average annual cost of the system, K, is the sum 
of (2-4) and (2-6). 

A, A I n CQ, A„A I 9 C[n-l]Q 1 

K = + -- V - + + CA + ^ . (2-7) 

Q-j^ 2 nQ^ 2 

The model has assumed that the unit price is independent 
of the quantity ordered. Therefore the C term in (2-7) is a 
constant, and can be dropped from (2-7) resulting in (2-8) 
being the average annual variable cost. 
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K 




+ 



I CQ 1 AX 

i + + 

2 + nQ 1 + 



I 2 C [n-l]Q 1 
2 



B. OPTIMAL POLICY 



( 2 - 8 ) 



Determining the optimal inventory policy involves finding 
the values of and n (call them Q* and n*) which minimize 
(2-8) for a given set of model parameters. To find this 
policy, rewrite (2-8) as 



K(Q 1 ,n) 






[I 1 + I 2 [n-1)]CQ 1 
2 



(2-9) 



Then since K(Q 17 n) is convex in Q 1 for fixed n, Q*(n) must 



satisfy the equation 



3K 

3Q, 



A , [I, + I,[n-l]]C 

t*i + Fr J r5 + ■ t - 0 ■ 

U 1 



or 



Q*(n) = 



2X [A, + -=•] 
1 n 



[I 1 + I 2 [n-1]] C 



1/2 



( 2 - 10 ) 



After substituting (2-10) for in (2-9) and collecting 
terms, (2-9) takes the following form: 



K(Q*(n),n) = {2X [A x + [I ] _ + I 2 [n-1] ] C} 



1/2 



or 



K(Q*(n),n) = 



2XA 1 I 1 C + 



2XA 2 I 1 C 

n 



+ 2X A, I 0 C [n-1] + 2X A_ I _C [— — — ] 
1 2 2 2 n 



1/2 



Note that the first term in (2-11) is equivalent to K 



w 



( 2 - 11 ) 



(Eq. (2-11) [Hadley and Whitin] ) . In fact, it is readily seen 
that each term in (2-11) contains a form of this classic 
formula. Let 



Kt . = 2X A • I . C 

ID 1 D 
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Then (2-11) can be rewritten 



K(n) = K (Q* , n) = {K 
1 w 



K 



21 



n 



K^ 2 [n-1] 



+ K 



22 






1/2 



( 2 - 12 ) 



It is immediately seen that (2-12) is positive for all 
feasible values of n (n >. 1) . Further, 



lim K (n) 
n-*l 




♦ . 



and for very large values of n, 
1. K (n) -*■ 



K 12 n 1/2 > 0 , 



2 . 



3. 



9K(n) 

9n 



9 2 K (n) 

. 2 

9n 



K 



12 



2n 

-K 



1/2 

12 



4n 



3/2 



> 0 , 



< 0 



If the assumption that n is an integer is temporarily re- 
laxed the slope of K(n) can be readily investigated by 
considering 



9K 

9n 



-K 



21 



K 



22 



n 



2 

n 

W 



2K (n) 



+ K 



2 

12 



(2-13) 



If (2-13) is positive when n = 1, then it follows that it 
is positive over all n, and the optimal value of n is n* = 1. 

If (2-13) is negative when n = 1, then there exist values of 
n such that K(n) is less than K(l). However, since it has 
been shown that (2-13) is positive for very large n, it follows 
that n* is finite, and satisfies 
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3K 

3n 



= 0 , 



or 



n 



* = 



K 2 - K 2 ‘ 
*21 *22 



K 



12 



1/2 



(2-14) 



A further consideration of the case where (2-13) is nega- 
tive at n = 1 reveals that 



3 2 K 

3n 2 



> 0 



over the range of n values from 1 to beyond n*. Let the value 
of n where 

& = 0 
3n 

A A 

be denoted by n. Then in the range 1 £ n £ n the function 

A 

K(n) is convex. For values of n > n, (2-12) is not convex. 
However, (2-13) will always be positive. Therefore n*, given 
by (2-14) , is the optimal non-integer value of n. The reader 
will note that the above arguments indicate that (2-12) is 
pseudo-convex . 

Since K(n) is pseudo-convex over all n, it follows that 
K(n) is also pseudo-convex on the integer values of n. If n* 
is not an integer, then the optimal integer value of n will be 
either n (the smallest value greater than n*) or n (the 
largest value less than n*) . The integer optimal will be the 
value of n corresponding to K* = min [K (n) , K (n) ] . 

In summary, the optimal inventory policy is found as 
follows : 
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1 . 



2 

Compute I< 2 ^/ 



If this inequality, holds then 




2 



Test K^, - K~ 0 1 K? 



n optimal =1. Go to step 8 



3. If - K ^2 < then compute n* using (2-14) 



4. If n* is an integer then n optimal = n*. Go to 
step 8. 

5. If n* is not an integer, compute n and n. 

6. Compute K(n) and K(n) from (2-12). 

7. K* = min {K(n , K (n) } and n optimal is the value of n 
corresponding to K* . 

8. Compute QJ from (2-10). 

C. NUMERICAL EXAMPLE 

As an illustration of the algorithm, consider the follow- 
ing problem where 



I 



1 



0.75 



A 2 = $200.00 



I 



2 



0.50 



X = 100 units/year 



A 1 = $25.00 C = $100/uni t 

2 

Begin by computing ^ = 2XA^I^C . 




Next compare - K^ 2 = 1.0 x 10 6 > K^ 2 = 2.5 x 10 5 
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so compute n* using (2-14) . The value is 2.0 which is 
integer. Therefore n optimal = 2. Using (2-10) the optimal 
order quantities are 

Q* = 166.67 units 

and 

Q* = 166.67 x 2 = 333.34 units. 

D. PARAMETER ANALYSIS 

The primary objective of this analysis was to observe the 
effect of parameter variations on optimal n. The values of 
four of the five parameters were fixed and the value of the 
remaining parameter was allowed to vary. The process was then 
repeated for each of the other parameters so that the effect 
of each could be observed. 

Representative graphs of optimal n versus a given para- 
meter are illustrated in Figures 7 through 10. In Figure 7 
it is seen that n optimal is inversely related to A-^. How- 
ever, as noted in Figure 8, there is a direct relationship 
between n optimal and A 2 * Further, a comparison of Figures 
7 and 8 reveals that, for the parameters selected, n optimal 
is much more sensitive to A 2 than to A^. 

Figures 9 and 10 indicate the relationship of n optimal 
to 1^ and I 2 , respectively. The results indicate that n 
optimal is directly related to 1^, inversely related to 
and that n optimal is more sensitive to a change in holding 
cost at echelon two than at echelon one. 
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An analysis of the effect of demand rate on n optimal 
revealed that n optimal was completely insensitive to X over 
the range 5 <. X <, 10,000. 
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III. FINITE PRODUCTION AT THE SECOND ECHELON 



An additional two echelon variation to the deterministic 
models treated [Hadley and Whitin] addresses the problem of 
finite production. This model would have application in any 
situation where the upper echelon manufactures as well as 
warehouses a given product. 



A. MODEL FORMULATION 

The model, shown graphically in Figure 3, assumes that 
echelon two produces the needed material at a constant con- 
tinuous rate of \p units per unit time (\p >. X ) . The cost A 2 is 
now considered to be a set up cost which is incurred each time 
a new production run is initiated. This set up cost is inde- 
pendent of the quantity produced. The rest of the assumptions 
are identical to those of the EOQ model. 

The model will, as in Section II, seek optimal values of 
and Q 2 which minimize the average annual system cost. This 
cost is the sum 

K = K 1 + K 2 , 

where is the average annual cost at the ith echelon. In 
Section II it was shown that 



K. 



A l^ 1 

qT“ + I 1 ! 00 ! * 



To develop K 2 it is necessary to define the cycle length 



as 



T = t. + T + T, , 
1 p d 



(3-1) 
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Echelon II 
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Figure 3. Two-Echelon Finite Production Model. 



where is the time from the beginning of the cycle until the 
start of production, T is the production length, and T, is 
the time required to deplete the inventory at echelon two 
after production stops (see Figure 3) . Note that 



T 

P 



and 



T 



d 



Ql _ Q x 
X \p ' 



x ' 



(C?2 (^2 — Qp) 

^ ^ - 



Define i as the number of times that the quantity 
demanded at echelon two during the production period, 
can be seen that 



Q-L is 

v It: 



t, + T 



(-) 



(3-2) 



where t = ^ — . 

As in Section II, the variable cost per cycle at echelon 
two, K 2 (C), will have the form 



t, + T + T, 
1 p d 

K 2 (C) = A 2 + I 2 C Ljj(t) dt 

fc l 



From Figure 4 it follows that 



(3-3) 




+ T 



P 



+ T 



d 



2 (t) dt = f (i) + g (i + 1) + h (i + 2) 



t 
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where f(i) is the area under the first i saw teeth, g(i + 1) 
is the shaded area shown in Figure 4 and h(i + 2) is the area 
remaining under the curve. 




Figure 4. Finite Production Model, Area Under the Inventory 
Curve at Echelon Two. 



From Figure 4 it can be shown that 



qI [i-1] [i-1] [i-2] 

f(l) = 2 $ + “ 2 + ^“2 



and 



(3-4) 



2X 



2X 



h(i) = 



Q x [n-i] [n-i-1] 

21 



(3-5) 
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From Figure 5 the area g(i + 1) is found to be 



g (i + 



Figure 5 



1 ) = 



-q\ [n-1] 2 Q^[i 2 ~l] 



2ip 



2X 



Q? [i [n-i] ] Q? { [i-1] [i^-X ] } 

4- i -f i 

X X2 



(3-6) 



!■<— a —4 




- lir-f 1 "- 1 " 

= — t»-l] - 4 [i-1] 

“ °1 [i-U[t^] 

= 

= Q 1 [n-1] - Q 1 f- [i-1] 

Finite Production Model. Area Under 
i + I s * 1 Saw Tooth. 
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Substituting (3-4) , (3-5) and (3-6) into (3-3) yields 



k 2 (c) 



I.CQ 2 I,CQ?if»[i-l] 

a j- £ i. + Z. x 

2 2l P 2 \ 2 



I 2 CQ 2 [n-i] [n-i-1] I 2 CQ 2 [n-l] 2 



2\ 



2X 



I 2 CQ 2 iMi-l] I 2 CQ 2 [i[n-i]] 

^ + X 

i,C[i-l] [ty-X ] 

+ — ^ . (3-7) 

X^ 

After dividing (3-7) by the total cycle length and collecting 
terms K 2 can be written 



K. 



A 2 X I 2 CQ 1 [n-l] I 2 CQiX [n-2] 
nQ^ + 2 2tJ) 



(3-8) 



The system average 'annual variable cost is found by summing 



and K 2 . 



A, [Ii + I? tn-1) - I 9 ^[n-2]]CQ, 

K(Q lf n) = [A l +~] + - — 2 ^ — - (3-9) 



B. OPTIMAL POLICY 



The development of an algorithm for determining the optimal 
inventory policy in the finite production model follows an 
argument that is identical to the one discussed in Section II. 
Since (3-9) is convex in Q 1 for fixed n, Q*(n) must satisfy 



9K (Q 1 , r.) 



0 , 



25 



or 



Q* (n) 



r A 2 

2 [A ,+— ]X 
In 


' 


C[I 1 +I 2 [n-l] - 


I 2 ?r [n - 2 1 



1/2 



When (3-10) is substituted for in (3-9) then 



K (n) = 



+K -yr + K i2 [[n ~ 1] " 



V 11 



K^ 2 [ [n-1] - ^[n-2]] 



n 



1/2 



(3 



(3 



where K7 . = 2XA.I.C. 

13 13 

It is apparent that (3-11) is positive for all values of n 
Further , 

lim K(n) = K 2 + K 2 + K 2 V + K 2 i. 

^ 11 21 12iJj 22 ty . 

and for very large values of n, 



2 . ^ 



3. 



K (n) - > -K 1 2 n 


9K _ 


K 12 


3n 


2n 1/2 


8 2 K 


_ - K 12 


9n 2 


An 3 '' 2 



> 0 , 



< 0 . 



If the assumption that n is an integer is temporarily 
relaxed the slope of K(n) can be investigated using 



BK 

3n 



2 2 
^21 2 X ^22 

-|i ♦ + -|i 



2XK 



22 



n 



2 K (n) 



n 

T7T 



\pn 4 



(3 



- 10 ) 



-ID 



- 12 ) 
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If (3-12) is positive when n = 1, it is also positive over 
all n and it follows that the optimal value of n is n* = 1. 
If (3-12) is negative when n = 1, then there exist values of 
n such that K(n) is less than K(l). However, since it has 
been shown that (3-12) is positive for very large n, n* must 
be finite and satisfy 



3K 

dn 



= 0 , 



or 



n 



* = 



2 ~2 
k 2 i - k 22 



n 1/2 



K 



12 



where 



~2 

K 22 



and 



K 



2 r ^-2A, 
22 l . i p J 



/ 



K 



2 

12 





(3-13) 



A further consideration of the case where (3-12) is nega- 
tive for n = 1 will show that 

> 0 
9n 

over the range of n values from 1 to beyond n*. Let the value 
of n where 




0 



be denoted by n, and it follows that in the range 1 <. n <. n 

K(n) is convex. For values of n > n K(n) is net convex, but 

♦ 
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(3-12) is positive. Therefore n* is the optimal non integer 
value of n and (3-11) is pseudo-convex. 

Since (3-11) is pseudo-convex over all n, it follows that 
it is also pseudo-convex on the integer values of n. So if 
n* is not an integer, then the optimal value of n will cor- 
respond to 

K* = min[K(n), K(n)] . 



The optimal policy algorithm can now be written. 

2 ~2 ~2 

1. Compute K 2 ^, K 22' an< ^ K 12 * 

2. If <. K^ 2 + K^ 2 then n* = 1. Go to step 8. 

3. If X-21 > K 22 + ^12 ^ en compute n* using (3-13). 

4. If n* is an integer then n optimal = n*. Go to step 



5. If n* is not an integer compute n and n. 



6. Compute K(n) and K(n) from (3-11). 



7. K* = min[K(n), K(n)] and n optimal is the value of n 
corresponding to K*. 

• * 

8. Compute Q* from (3-10). 

C. SPECIAL CASES OF \p 
When ip = X 



K (n) 



Q x 2 nQ^ 2 ' 



and 
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(3-14) 



lim K(n) 
n+°o 



°r + — 5 — 



Thus it is seen that if the production rate is equal to the 

demand rate, the optimal policy is to start the production 

line and never let it stop. 

As one would expect, when ip -*■ 00 

A X I -I CQ. A X 

lim K(n) = 5 ^ + - T - + 

I 2 CQ 1 [n-l] 

+ 2 ' 

and it is immediately seen that this model degenerates to 
the EOQ model of Section II. 



D. NUMERICAL EXAMPLE 

Consider the problem where 



I 1 = 0.50 



I 2 = 0.20 



X = 25 units/year 

ip = 75 units/year 



C = 

A 1 = 
A„ = 



$100/unit 
$50/uni t 
$200/unit 



First compute Kf ^ , 


K 22 and Kj 2 . 










K 2 - 

K 21 


5.0 x 


10 5 


K h " 


1.25 


X 


10 5 


K 2 - 

K 22 


2.0 x 


10 5 


*=22 = 


6.67 


X 


o 

i — i 


2 




4 


^2 






4 


K 12 = 


5.0 x 


10^ 


K 12 “ 


3.33 


X 


10 
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Next compare with — K 21 = x > ^22 + ^12 = x 

so compute n* using (3-13). The value is 3.64 which is not an 
integer. Therefore n = 4 and n = 3. 

From (3-11) 

K (n) = $685.00, 

K (n) = $694.00, 

and n optimal is 4. Using (3-10) 

Q* = 7.2 units 

and 

Q* = 28.8 units 

E. PARAMETER ANALYSIS 

Figures 11 and 12 illustrate the effect of production rate 
on optimal n. If the holding cost at echelon two is small 
then optimal n was relatively insensitive to a change in 
production rate (Figure 11). However, the results of Figure 
12 show that this sensitivity increases as increases. 

The responses of this model to parameter variations were 
identical to the responses observed in Section II. 
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IV. THE TWO ECHELON MODEL WITH BACKORDERS PERMITTED 



A. MODEL FORMULATION 

The EOQ and finite production models were based on the 

assumption that all demands would be immediately satisfied. 

This assumption will now be relaxed in order to investigate 

a system in which all demands are ultimately satisfied, but 

where it is permissible to accumulate backorders at the lower 

echelon. No backorders are allowed in echelon two. 

In this model demands which occur when the lower echelon 

is out of stock are backordered against future procurement. 

When the procurement arrives these backordered demands are 

met, and any excess quantity is placed in stock. Each back- 

/\ 

order at echelon one results in a cost of the form it + ir, 

/s 

where i is a fixed change and it is a time dependent cost. 

As noted, if there were no costs associated with incurring 
backorders [Hadley and Whitin] , then it would be optimal to 
never have any inventory on hand. Conversely, if the back- 
orders are sufficiently expensive, then the optimal policy 
would be to never incur any. However, for an intermediate 
range of backorder costs, it will usually be optimal to incur 
some backorders toward the end of the cycle. 

With the exceptions noted in the preceding paragraphs, the 
model is predicated on the same assumptions used in Section II. 

Define s as the number of backorders at echelon one when 
an order is received. Then from Figure 6 it follows that the 
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on hand inventory at echelon one immediately after receipt of 
an order is (Q^ - s) units. It is assumed that >. s. Then 

the echelon one on hand inventory varies from (Q^ - s) to 
zero. 



Echelon II 
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Figure 6. Two-Echelon Backorders Allowed Model. 



The optimal policy continues to be that which minimizes 
average annual variable system cost subject to the model 
constraints. Since the model does not permit backorders at 
the second echelon, the average annual variable cost at 
echelon two will be identical to equation (2-4), which for 
convenience is restated here. 
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K, 



A 2 X I 2 CQ 1 [n-lJ 
Kq7 + 2 



(4-1) 



Development of the echelon one variable costs follow 
identically the single echelon backorders permitted model of 
Hadley and Whitin. In the interest of brevity, this develop- 
ment will be omitted and the results (equation 2-17) [Hadley 
and Whitin] merely stated. 



is 



A X ijCK^-s] ; s 2 

K 1 - + ' Sq-" - + q w.fi- 



(4-2) 



Summing (4-1) and (4-2), the average annual variable cost 



A i A V 1 i c [ Q r s ]' 

K = Q]“ + + 2Q^ 



, 1 r 'i , TTS , 

+ [1TXS + — 1 1 



(4-3) 



I^CQ^n-l] 



B. OPTIMAL POLICY 

The optimal solution seeks the values of the decision 
variables Q*, s*, and n* which minimize (4-3) subject to the 
constraint that all demands are ultimately satisfied. To do 
this, note that (4-3) can be rewritten as 



where 



K = 



X = 



AX 



IC ^1 T r .ttXs Xs^ 

~ - i i cs + or + 207- ' 



I^C + TT , 



(4-4) 



and 



A = 



A, + — 
1 n 



33 



I = l 1 + l 2 [n-'l] 



Then from 



M. = 0 

as ' 



it follows that 



s = 



X l CQ l 



- ttA 



X 



(4-5) 



or 



ttA + Xs 



J 1 C 



Next, 

8K 

3Q X 

results in 
Qi 



-2AA 


IC 2ttAs 


„ 2 + 


2 ^ 2 


2Q X 


2Q 1 


r * 


2 


2AA + 


2ttAs + Xs 




IC 



Xs 

2Q! 



= 0 



ni/2 



(4-6) 



(4-7) 



To solve for Q£ and s* in terms of n and the system parameters, 
substitute (4-7) into (4-5) and (4-6). After collecting terms, 
the resulting equations are 

11/2 



Q*(n) = 



2 2 

2A X - it X 
2 2 

ICX - C 



(4-8) 



and 



s*(n) = 



X 



2AA X - tt 2 A 2 



2 2 

ICX - i x C 



1/2 



ttA 

X 



(4-9) 



The expressions (4-8) and (4-9) can now be back substituted 
for and s in (4-4) to yield an expression for the average 
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annual variable cost, K, as a function of n and the system 
parameters. After regrouping the terms this expression can 
be written 



K(n) = 



1/2 ^ 1/2 

[IC [I^C+tt] — I 2 C 2 ] [2AX[I 1 C+7T]-tt 2 X 2 ] I-jCttX 

/v T P+ A 

x l'- + 7r 



[I-jC + it] 



(4-10) 



Equations (4-8) , (4-9) , and (4-10) are valid only if 

2 2 

2AXIC - it X >.0. Otherwise an analysis similar to that of 



Hadley and Whitin is required for the case when tt and tt are 
positive. Such an analysis is beyond the scope of this thesis. 

Because of the complex form of K(n) given by (4-10) as a 
function of n no attempt was made to evaluate whether or not 

A 

K(n) was convex or pseudo-convex in n when tt and tt are both 
positive. However, an example is presented later in which an 

A 

optimal value of n > 1 is obtained for tt and tt positive. 

If tt = 0 (4-10) reduces to 



I 1 C 

K (n) = [2AXIC - 2AXI 1 C 



1/2 



(4-11) 



To expand (4-11) , let 



K7 . = 2X a . I .C 
1 3 i D 



(4-12) 



Then , 



12 2 ^ 21 ^ 
K,n) = ^7? K ii x + V”- 1 ' + 



K 22 X[n 1) >2 T ^ K 21 I 1 c 

n 11 1 n 
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and since x = 



I^C + tc, K(n) reduces to 



K(n) = 



2 2 

^11 77 ^21 71 2 
- ~ T + — T + K‘ 2 [n-1] 

I C+u n[I C+i] 

i 



K 2 2 [n-1] 



n 



1/2 



(4-13) 



Following the analysis of Section II, it is immediately 
seen that (4-13) is positive for all feasible values of n, and 
that 



lim K (n) = 

n->l 



2 ~ 2 ~ 
*<11” + ^l 77 

I^C + TT 



1/2 



Also, as n gets very large 

1. K(n)+K 12 n 1/2 > 0 , 



2. 



9K . K 12 



8n 2n l/2 



> 0 , 



9 2 K (n) ^ K 12 „ A 

3 ‘ ^72 • 



Assume, for the sake of argument, that n is continuous 
for all n > 1. Then 



-K 



21 



9 2 K 

9n 2 



[I^C+tt] n 



2 K 22 
_l i<f r f f 

2 *12 2 



n 



K (n) 



1 / 2 " 



(4-14) 



and it is obvious that for n = 1 (4-14) is non negative if and 
only if 
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(4-15) 



I C 

K 21 - K 12 + K 22^ + ~ ^ ^ * 

IT 

If (4-15) holds, then (4-14) is positive for all n. This 
implies that over all n, the slope of K(n) is positive, and 
n* = 1. 



IF (4-14) is negative at n = 1, then there exist values 
of n such that K(n) < K(l). However, it has been shown that 
for very large n, (4-14) is positive, and it follows that n* 
is finite and satisfies 



3K 

3n 



0 , 



or 



n 



* — 



K 21 tt K 2 2 ^1 C + ^ 



1/2 



K 12 [I 1 C + v] 



(4-16) 



To complete the argument, note that when (4-14) is 
negative at n = 1, 




0 



over the range of n values from 1 to beyond n*. Let the value 
of n where 




3n 2 



be noted by n, and it follows that K(n) is convex over the 

✓S A 

range 1 <_ n <. n. For values of n > n, (4-13) is not convex, 
but in this range (4-14) is always positive. Therefore n* 
given by (4-16) is the optimal non integer value of n. The 
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reader should again note that the above arguments indicate 
that (4-13) is pseudo-convex. 

Since K(n) is pseudo-convex over all n, it is also pseudo- 
convex on the integer values of n. Therefore, if (4-16) does 
not yield an integer value then the optimal value of n can be 
found by evaluating K(iT) and K(n) as in Section II. 

An algorithm for finding the optimal policy (in the case 
where tt = 0) can now be stated. 

2 2 2 

1. Compute K^, K 22' anc ^ K 12 * 



(4-16) . 

4. If n* is an integer then n optimal = n*. Go to step 8. 

5. If n* is not an integer compute n and n. 

6. Compute K(n) and K(n) from (4-13). 

7. K* = min[K(n) and K(n)] and n optimal is the value of 
n corresponding to K*. 

8. Compute Q* from 




] . If the inequality 



TT 

holds then n optimal is 1. Go to step 8. 




2A [A n + — ] 
1 n 



1/2 




C[I 1 I 2 C[n-l] + I 1 tt + I 2 n[n-1]] 



and s* from 



s* (n) 



I 1 C Q* 1/2 



I^C + TT 
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C. NUMERICAL EXAMPLE 



Consider the problem where 



J 1 - 


0.75 


7T = 


$5.00 


I 2 


0.25 


A 

TT = 


$100 . 00/year 


A 1 = 


$25.00 


X = 


100 units/year 


A 2 " 


$100.00 


c = 


$100/unit 


Begin 


2 

by computing 


= 2XA.I 

l 


.C. 

3 



= 1. 5 x 10 6 

K^ 2 = 5.0 x 10 5 

K^ 2 = 1.25 x 10 5 

I c 

2 2 2 1 

Next compare K 2 ^ < [K^ 2 + K 22 ] [1 + ] . The inequality does 

ir 

not hold so compute n* using (4-17). The value is 1.69 which 
is not an integer. From (4-13) 

K (n) = $735.00 

and 



K (n) = $620.00 

Therefore optimal n = 1. The values of Q£ , Q* , and s* 
are now found to be 



Q| = 12.9 units 

= 12.9 units 

s* = 8.16 units. 
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* 






D. PARAMETER ANALYSIS 

As noted previously, computational difficulties preclude 
an analytical treatment of the general form of the. backorder 
model. However, the response of the model to parameter 
variations can be approximated. This can be accomplished by 
considering separately the two functions (3-1) and (3-2) and 
ignoring, for now, the integer restriction on n. 

Equation (4-1) can be rewritten 



K, 



V . I 2 C|Q 2 - Ql ] 

rr ~r 



and since the second term is never negative, it is obvious 

Using the Hessian it can be shown 
that (4-2) is convex in and s if and only if 



that K 2 is convex in Q 2 



Then 



2A 1 X [I^C + it] >. [ttX ] 



K 



AjA 



I 1 C[Q 1 - s] 



+ — + ^[ttXs + --] 



TTS 



(4-18) 



+ A 2^ + I 2 C ^ Q 2 ~ Q l^ 



Q. 



(4-19) 



is convex if and only if (4-18) holds. 

It can be shown that when (4-19) is not convex the optimal 
inventory policy will be realized when backorders are not 
allowed (s* = 0). Thus if (4-18) does not hold, the problem 
can be solved by setting s* = 0 and using the algorithm 
described in Section II to get Q| and n*. 
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For the case where (4-19) is convex, optimal Q^, Q 2 and 
s can be found from 



3K 

3Q. 



3K 

3Q. 



3K 

3s 



= 0 



(4-20) 



Equation (4-20) will be satisfied when 



Q* 



2X [nA^ + A 2 ][I^C + tt] - niT 2 X 2 
ntljC + tt + I 1 I 2 C 2 [n-l] + I 2 CTT[n-l] 



1/2 



(4-21) 



s* (n) = 



I-jCQjfn) - trX 

I^C + TT 



(4-22) 



and 



Q*(n) 



= nQ^(n) , 



(4-23) 



and the optimal policy can be found from the following 
algorithm. 

1. Set n = 1 . 



2. Compute Q*(n), s*(n), Q* (n) and K(n) using equations 
(4-21) through (4-23) and (4-19). 

3. Set n = 2 and repeat step 2. 

4. If K ( 1) < K(2), stop; n* = 1 and the optimal policy 
is known. 



5. If K(2) < K (1) , set n = 3 and repeat step 2. 

6. Continue solving K (n) by increasing n in steps of 1 
until K(n+1) >. K(n). Stop as soon as K(n) < K(n+1); n* = 
current value of n and the optimal solution has been found. 



41 



The algorithm described in the preceding paragraph was 
used to observe the effect of parameter variations on optimal 
n. In general, the results were similar to those observed in 
the EOQ model. For example, n optimal varied inversely with 
A^ and I 2 , but directly with A 2 and I 1 . 

The sensitivity of n to the various parameters was highly 
dependent on the selected values of the backorder costs. For 
example, Figure 13 illustrates the relationship between optimal 
n and I 2 in a situation where the backorder cost is small. 

Figure 14 illustrates this same example with a high time de- 
pendent backorder cost. The difference in the sensitivity is 
obvious . 

Figures 15 and 16 illustrate the relationship of optimal n 
to a change in A^. In Figure 15 the backorder costs were low, 
and optimal n was insensitive to a change in A^. In this 
example optimal n remained at two, and s* increased as A 2 
increased. When the time dependent backorder cost was high 
optimal n varied directly with A 2 (Figure 16) . However, the 
change in optimal n was not nearly as large as in the EOQ model. 

Figures 17 through 20 illustrate the relative magnitude of 
Q* , Q* , and s* as A 2 and I 2 are allowed to vary. The breaks 
in these curves occur at points where optimal n changes. 
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V. RECOMMENDATIONS FOR FURTHER STUDY 



The analysis of the backorders permitted model should be 

A 

completed for all values of the backorder costs it and tt in- 

/A 

eluding the case where tt > 0 and tt = 0. Hadley and Whitin 

discuss this case in their single echelon development and 

/\ 

conclude that if it = 0 then s* is either 0 or infinite. 

Initial investigations indicated that in the two echelon model 

-A 

there could exist a finite value of s* > 0 when tt = 0. How- 

ever, this investigation was not completed and no conclusions 
were reached. 

It is felt that the analytical argument used throughout 
the paper could be applied to the general backorder case where 
tt and tt are both greater than zero. However, the complexity of 
the equations would greatly complicate this development. 

The models discussed in this thesis should be extended to 
more than one activity at each echelon. If it can be assumed 
that all activities within a given echelon order at the same 
time, then it is particularly easy to include multiple activ- 
ities within an echelon. For example, assume that there are 
K activities in echelon one. Then in the EOQ model would 
take the form 



K. 



k 

l 

i=l 



A, . A . I. .CO, • 
li l li li 



li 



If the values of the parameters A^ , A^, 1^, and C are constant 
for all i, i = l,...,k, then the expression for reduces to 
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the trivial case of 



K' = k K ± , 

where is found in Section II. The expression for 
remains unchanged. 

The problem of multiple activities at the top echelon is 
usually not a very interesting one, since each top echelon 
activity is usually responsible for supplying specified activ- 
ities at the next lower echelon. Therefore, unless the model 
provides for lateral resupply actions, the problem of multiple 
activities at the top echelon can be reduced to the sum of a 
series of independent problems, each containing one top 
echelon activity. 

A related problem which should be investigated is the ex- 
tension of the models to more than 2 echelons. In the case 
where each echelon is limited to one activity, this extension 
would not be difficult. Indeed, given the assumption that there 
is an integer ratio between the quantity ordered by successive 
echelons, the model would probably take the form 
m 

K - I K • , 

3=1 D 



where the relationship between the Q^'s is 



= n. 



= n. 



Q 2 n 2 



n. 



Qj n j_ x * Qj_i , n , n i Q i ' 



j-i 

n 

i=l 



and nK is an integer for all k, k = 1 



, j-1. Because of the 
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integer property of the 1 s , the solution technique outlined 
in Chapter II could be used to solve recursively for K. 

i i_ 

If the m-echelon problem allows activities at the 1 
echelon, i = l,...,m, then for large values of and m, the 
sheer magnitude of the problem would make its solution ex- 
tremely difficult. In fact, it is doubtful that the proce- 
dures suggested in this paper could be utilized for a problem 
of this type. However, if the model is small enough, one 
could at least get a "feel" of the model's behavior. 
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VI . SUMMARY 



Three of the deterministic models of Hadley and Whitin 
have been restructured as two echelon models. Equations for 
the average annual variable cost are derived, and standard 
mathematical programming techniques are utilized to find the 
optimal inventory policy. The optimal policy is defined as 
that which minimizes average annual variable cost subject to 
the constraints of the model. 

General solutions are obtained for the EOQ and finite 
production models. Because of the complexity of the cost equa- 
tion the analytical solution of the backorders permitted model 
is limited to the case where it = 0 . However, a technique is 
developed which can be used to find the optimal policy in the 
general case. 

A parametric analysis is conducted for each model in order 
to study the behavior of the models under parameter variations. 
The behavior of the curves of n optimal for the EOQ and finite 
production models were identical except that the finite produc- 
tion model enjoyed a significant cost advantage. The behavior 
of the backorder model was very dependent on the magnitude of 
the backorder cost. For example, if the model had a low time 
dependent backorder cost then parameter variations had little 
effect on the optimal behavior of the model. However, if the 
backorder cost was set sufficiently high then the model's 
behavior v/as similar to that of the EOQ model. Figures 7 
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through 20 illustrate the behavior of the three models under 
varying parameter values. 
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Parameter Values 
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Figure 7. Optimal n : A, • EOQ Model. 
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Figure 8. Optimal n : A- • EOQ Model. 
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Figure 9. Optimal n : I, • EOQ Model. 
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Figure 10. Optimal n : I~ • EOQ Model. 
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Figure 11. Optimal n : Production Rate. Finite Production Model 
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Figure 12. Optimal n : Production Rate. Finite Production Model 
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Figure 13. Optimal n : I 2 with Low Backorder Cost. Backorders Allowed Model 
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Figure 14. Optimal n : I2 with High Backorder Cost. Backorders Allowed Model. 
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Figure 15. Optimal n : A_ with Low Backorder Cost. Backorder Allowed Model 
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Figure 16. Optimal n : A 2 with High Backorder Cost. Backorder Allowed Model. 
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Q* , and s* : Backorder Model 

with Low Backorder Costs. 
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Figure 18. Q* , Q* and s* : Backorder Model v-;ith 

High Backorder Costs. 
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and s* : Backorder Model with High Backorder Costs. 
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